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Abstract: Due to their special phylogenetic position in the Euarchontoglires and close affinity to primates, tree 
shrews have been proposed as an alternative experimental animal to primates in biomedical research. However, the 
population genetic structure of tree shrews has largely remained unknown and this has hindered the development of tree 
shrew breeding and selection. Here we sampled 80 Chinese tree shrews (Tupaia belangeri chinensis) in Kunming, China, 
and analyzed partial mtDNA control region sequence variation. Based on our samples and two published sequences from 
northern tree shrews (T. belangeri), we identified 29 substitutions in the mtDNA control region fragment (~ 604 bp) 
across 82 individuals and defined 13 haplotypes. Seventeen samples were selected for sequencing of the cytochrome b 
(Cyt b; 1134 bp) gene based on control region sequence variation and were analyzed in combination with 34 published 
sequences to solidify the phylogenetic pattern obtained from control region data. Overall, tree shrews from Kunming have 
high genetic diversity and present a remarkable long genetic distance to the two reported northern tree shrews outside China. 
Our results provide some caution when using tree shrews to establish animal models because of this apparent genetic 
difference. In addition, the high genetic diversity of Chinese tree shrews inhabiting Kunming suggests that systematic 
genetic investigations should be conducted before establishing an inbred strain for medical and biological research. 
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Tree shrews (Tupaia belangeri) are small mammals 
of a squirrel-like appearance formerly placed in the 
Primates order despite a lack of derived features 
characteristic of primate species (Sargis, 2004). It is 
currently placed in the family of Tupaiidae, Scandentia, 
which together with Dermoptera and Primates, form the 
group Euarchonta (Murphy et al, 2001; Helgen, 2005). 
Although the exact phylogenetic position of tree shrews 
remains resolved, this animal is still regarded as a close 
relative to primates (Nie et al, 2008; Peng et al, 1991; 
Seiffert et al, 2003). Given that tree shrews are easy to 
raise, cheap to maintain, have a high reproductive rate, 
small body size, and close affinity to primates, this 
animal may be an alternative to the use of primates in 
biomedical research. For example, tree shrews have been 
used in creating animal models for hepatitis B virus and 
hepatitis C virus infection (Amako et al, 2010; Cao et al, 
2003; Köck et al, 2001; Ren & Nassal, 2001; Su, 1987; 
Xu et al, 2007; Yan et al, 1996). Tree shrews were also 
successfully used in the study of pathogenesis of human 
myopia and psychosocial stress diseases (Cao et al, 2003; 
Fuchs, 2005; Norton et al, 2006). In addition, tree shrews 
have the highest brain to body ratio of all mammals, 
including humans (Peng et al, 1991; Previc, 2009), which 
raises the possibility of investigating brain function and 
creating animal models for human neuro-degenerative 
diseases. 

There are four genera within the family Tupaiidae, 
consisting of Tupaia, Anathana, Urogale, and 
Dendrogale (Helgen, 2005). In genus Tupaia, a total of 
15 species are recognized and broadly distributed across 
Southeast Asia 
Philippines, Java, Borneo, Sumatra and Bali (Olson et al, 
2005; Peng et al, 1991). Chinese tree shrews (Tupaia 
belangeri) are distributed across southwest China. 


including southern China, India, 


According to geographical distribution and 
morphological characteristics, Wang (1987) proposed 
that Chinese tree shrews could be divided into six 
subspecies: Tupaia belangeri chinensis distributed in 
most of parts of Yunnan, Yun-Gui plateau, and 
southwest Sichuan; T. b. gaoligongensis distributed in 
central and north Gaoligong Mountain; T. b. modesta on 
Hainan Island; T. b. tonquinia distributed in southwest 
Guangxi; T. b. yaoshanensis distributed in northwest 


Guangxi; and T. b. yunalis distributed in central and 
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south Yunnan, northwest Guangxi and southwest 
Guizhou. Recently, Jia et al (2009) studied the 
morphometrics of the skull and mandible of T. b. 
chinensis collected in Luquan, Jianchuan, Lijiang, and 
Yunlong Counties of Yunnan and morphological 
differences between localities. In combination these 
studies suggest that Chinese tree shrews show regional 
differentiation. 

Despite growing attention and the strong potential 
of tree shrews as an animal model for human disease, 
scant biological data about these species have been 
collected, especially at the molecular level (Peng et al, 
1991). This deficit has constrained the utilization of tree 
shrews as an animal model in biomedical studies. 
Although the critical problem regarding laboratory 
breeding and domestication of wild tree shrews was 
solved more than two decades ago (Peng et al, 1991), an 
inbred strain of tree shrews has not been established 
because, among other factors, of a lack of genetic data. 

Here we aimed to analyze the genetic structure of 
tree shrews collected from urban Kunming, Yunnan, 
China to provide useful information to establish several 
inbreeding lines from these captured animals. We used 
mitochondrial DNA (mtDNA) as a genetic marker as it is 
useful in revealing the origin and demographic history 
and overall genetic diversity of domestic and wild 
species (He et al, 2008). 


1 Materials and Methods 


1.1 Sampling, DNA amplification and sequencing 

We collected ear or muscle tissue from 80 Chinese 
tree shrews (T. b. chinensis) inhabiting the suburb of 
Kunming, Yunnan, China. Animals were raised at the 
experimental animal core facility of the Kunming 
Institute of Zoology and procedures were approved by 
the Institutional Animal Care and Use Committee of the 
Kunming Institute of Zoology, the Chinese Academy of 
Sciences. 

Genomic DNA was extracted from tissue samples 
using the standard phenol/chloroform method. A 
fragment with a size of 1488 bp that covers the entire 
mtDNA control region of Chinese tree shrews was 
amplified using primer pairs L15356: 5'-CTCAAGGAA 
GAAGAACAATA-3' / H50: 5'-TGTATGTTTATGGAG 
TCTATG-3'. Primers were named relative to their 
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positions in the complete mtDNA sequence of northern 
tree shrew (GenBank accession number NC_002521). 
PCR amplification was performed in a volume of 50 uL 
reaction mixture containing 100 ng of DNA, 10 mmol/L 
Tris-HCl (pH 8.3), 2.5 mmol/L MgCl, 50 mmol/L KCl, 
10 pmol/L of each primer, and one unit of LA Taq 
polymerase (Takara, Dalian, China) following 35 cycles 
of 50 s at 94°C, 50 s at 55°C and 60 s at 72°C. The 
complete mitochondrial Cyt b gene was amplified and 
sequenced in 17 Chinese tree shrews representing 
different haplotypes as identified by mtDNA control 
region data. We used primer pair L14058: 5'- 
ACAAAAACACAATATAAGGCA -3' / H15298: 5'- 
CTCCGTTTCTGGTTTACAAG -3' to amplify the Cyt b 
gene, applying the same amplification conditions and 
procedures as used for mtDNA control region. PCR 
products were purified on spin columns and directly 
sequenced using Big Dye Terminator v3.1 Cycle 
Sequencing Kit (Applied Biosystems, Carlsbad, USA) on 
an ABI 3730 DNA sequencer according to the 
manufacturer's manual. Two inner sequencing primers 
(H16093: 5'-CAAGTTAAGTCCAGCTACAAT-3', for 
mtDNA control region; L15059: 5'-TCCCATTCCTT 
CACACGTC-3', for Cyt b gene) and primers L15356, 
L14058, 
sequencing of the analyzed region. The mtDNA 


and H15298 were used for overlapping 


sequences of Chinese tree shrew have been deposited in 
GenBank with accession numbers HQ836245- HQ836324 
(for mtDNA control region sequences) and HQ836325- 
HQ836341 (for Cyt b gene sequences). 

1.2 Data Analysis 
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We analyzed a fragment of mtDNA control region 
(~604 bp) in 80 Chinese tree shrews, corresponding to 
15416-16017" region in mitochondrial genome of T. 
belangeri (GenBank accession number NC_ 002521). 
Two reported individuals of northern tree shrews (Tupaia 
belangeri) outside China were retrieved from GenBank 
and were aligned with those of Chinese tree shrews (Tab. 
1 and Fig. 1). Sequence variations were exported using 
Mega 4.0 (Tamura et al, 2007). The available Cyt b gene 
sequences of tree shrews in GenBank, including 32 
Chinese tree shrews and two northern tree shrews outside 
China (Tab. 1), were also analyzed together with 17 
newly sequenced samples. Relationships among the 
lineages were explored using the median-joining 
networks (Bandelt et al, 1999). 
constructed manually and were checked using Network 


Networks were 
4.1 (http://www. fluxus-engineering.com/sharenet.htm), 
in which transitions and transversions were equally 
weighted. Nucleotide diversity (a) and haplotype 
diversity (A) were estimated according to Nei (1987) and 
using DnaSP 5.0 (Librado & Rozas, 2009). 

To further determine the phylogenetic position of 
Chinese tree shrews we retrieved the Cyt b gene 
sequences of Tupaia glis (AY321645), T. longipes 
(AY321651), T. salatana (AY321654), and Dendrogale 
melanura (AY321634, as the outgroup) from GenBank 
and constructed a neighbor-joining (NJ) tree using Mega 
4.0 (Tamura et al, 2007). The maximum composite 
likelihood nucleotide substitution model was employed. 
The rate variation among sites, with four rate categories 
for the nucleotide gamma shape parameter, was 


Tab. 1 mtDNA sequence information for tree shrews analyzed in this study 














Species GenBank accession number No. h° T°’ References 
Control region 
Chinese tree shrews HQ836245- HQ836324 80 0.790 + 0.032 0.00436 + 0.00017 This study 
Northern tree shrews outside AF217811" 1 _ _ Schmitz et al, 2000 
China AJ421453 1 = = 
Cyt b gene 





HQ836325- HQ836341 


Chinese tree shrews AY321635-AY321637 


EU376531, EU531749-EU531776 


b 
Northern tree shrews outside AF217811 


China AJ421453 








17 0.926 + 0.035 0.00338 + 0.00029 This study 
3 From GenBank 
29 Jia et al, 2008 


1 of 2: 


1 = y 


“Nucleotide diversity (x) and haplotype diversity (h) were estimated according to Nei (1987) . We pooled all Cyt b sequences of Chinese tree shrews together to 


estimate these two parameters. 


> This sequence has been curated by NCBI staff and has the accession number NC_002521. 
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Fig. 1 Sequence variation of mtDNA control region fragments 
(~604 bp) of 80 Chinese tree shrews and two reported 
northern tree shrews outside China 

The number of individuals sharing the same haplotype are listed in the right 

column, under N. Dots (-) denotes identical sites. The numbering of each 

position was based on the mitochondrial genome of the tree shrew 

(GenBank accession number NC_002521). 


estimated using jModelTest 0.1.1 (Guindon & Gascuel, 
2003) which finally yielded the alpha value 0.64. 


2 Results 


2.1 Sequence variation and haplotype distribution 

The mtDNA control region sequences of tree 
shrews presented length mutations, which were caused 
by insertion/deletion (indels) of tandem repeat or length 
mutations of a stretch of duplicated nucleotide in region 
15449-15456. One sample (TS07) had an extra insertion 
of CATT repeat in region 15557-15564 which contained 
two tandem repeats in the analyzed fragment (~604 bp). 
After the indels were excluded, we found 29 sequence 
variations in the analyzed fragment of mtDNA control 
region, which defined 13 haplotypes (HI—H13; Fig. 1). 
Among them, haplotype H5 was dominant and occurred 
in 30 individuals, followed by haplotypes H1 (in 14 
individuals), H2 (in 12 individuals) and H12 (in 8 
individuals). The other haplotypes had a low frequency 
and were shared by no more than four samples. The two 
reported northern tree shrews outside China shared 
haplotype H13, and presented strong differences to the 
haplotypes identified in Chinese tree shrews (Fig. 1). The 
overall genetic diversity of 80 Chinese tree shrews (x = 
0.00436 + 0.00017 and h = 0.790+0.032) was relatively 
high. 

One hundred and twelve variable sites were found 








in 51 Cyt b gene sequences (including 49 Chinese tree 
shrews and the two northern tree shrews outside China), 
which defined 26 haplotypes (hl—h26). Among these 
sequence variations, 17 substitutions caused amino acids 
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changes. Some representative samples bearing different 
mtDNA control region haplotypes, e.g. HS and H7, were 
found to share identical Cyt b sequence of h4. However, 
the two randomly selected individuals sharing mtDNA 
control region haplotype H5 differed from each other by 
one substitution in the Cyt b sequences (h4 and h5). 
Samples bearing unique mtDNA control region 
haplotypes (H1, H2, H3, H8, H9, H12, and H13) also 
differed from 
substitutions in the Cyt b gene. Remarkably, the two 


each other by some nucleotide 
northern tree shrews outside China shared the same Cyt b 
sequence and presented at least 89 mutation distances to 
Chinese tree shrews (Fig. 2). The nucleotide diversity of 
49 Chinese tree shrew Cyt b gene sequences was 
0.00333 + 0.00032. 
2.2 Network profiles of mtDNA haplotypes in 
Chinese tree shrews 
Fig. 3a shows the network file of mtDNA control 
region haplotypes in 82 tree shrews analyzed in this 





study. The network of all Chinese tree shrews followed a 
roughly star-like pattern, with haplotype H3 residing in 
the central location and being the nodal type. However, 
the frequency (5%, 4/80) of H3 was quite low in our 
samples. All other haplotypes differed from H3 by no 
more than a 3-mutation distance, suggesting relatively 
recent divergence. The three dominant haplotypes (H1, 
H2, and H5) were distributed at the terminal branches of 
the network. Up to 18-mutation distance existed between 
Chinese tree shrews (haplotypes H1-H12) and the 
reported northern tree shrews outside China (haplotype 
H13). 

The network profile of the Cyt b haplotypes 
followed a similar pattern observed for mtDNA control 
region sequences (Fig. 3b). Because we selected the 
samples for Cyt b sequencing based on the haplotypes 
identified by control region sequence information, an 
estimation of the frequency of each Cyt b haplotype may 
be biased, despite a fact that we took 32 published 
Chinese tree shrews samples into consideration. 
Nonetheless, the phylogenetic relationship as indicated 
by the network of Cyt b haplotypes might offer further 
information for us to discern the genetic difference 
between different lineages. Consistent with the pattern of 
mtDNA control region sequence data, haplotype h3 was 
also located in the central location of the network, 
suggesting its ancestral status. The difference of Cyt b 
gene nucleotide changes between Chinese tree shrews 
and the reported northern tree shrews (h26) was 
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Fig. 2 Sequence variation of mtDNA cytochrome b gene fragments (1134 bp) of 17 Chinese tree shrews newly 


sequenced in this study and 34 individuals from GenBank 
The detailed sources of these reported sequences are listed in Tab. 1. Dots (-) denotes identical sites. Substitutions were numbered according to the 


mitochondrial genome of northern tree shrew (GenBank accession number NC_002521). The number of individuals sharing certain mtDNA control region 


haplotypes (in parentheses) is listed in the right column. 
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Fig. 3 Network profiles of the mtDNA control region (a) and Cyt b gene (b) haplotypes of tree shrews 
The links are labeled by nucleotide positions to designate variations. The order of the mutations on a branch is arbitrary. The area of each circle is proportional 


to haplotype frequency. White color denotes individuals analyzed for both mtDNA control region and Cyt b gene sequence variation. Chinese tree shrew Cyt b 


gene sequences from GenBank are marked black. 
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remarkable, with a ratio (89/1134) even higher than that 
of the mtDNA control region (18/604). 
2.3 Phylogenetic position of Chinese tree shrews 

We constructed an NJ tree of the 26 Cyt b 
haplotypes of Chinese tree shrews, together with four 
species from GenBank (Tupaia glis, Tupaia longipes, 
Tupaia salatana, and Dendrogale melanura), to further 
demonstrate the affinity of Chinese tree shrews to other 
species (Fig. 4). All Chinese tree shrews haplotypes were 
clustered together with high bootstrap support. Chinese 
tree shrews are clearly separated from northern tree 
shrews outside China and other species, but presented an 
overall closer relationship to northern tree shrews. 
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h23 
h1 

79) h24 
h22 


M 19 Chinese tree shrews 


h18 (Tupaia belangeri 
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Fig. 4 The neighbor-joining (NJ) tree of the Cyt b gene of 26 
Tupaia belangeri haplotypes and four other species (T. 
glis, T. longipes, T. salatana, and Dendrogale 
melanura) 

Tree was rooted by Dendrogale melanura. The values on the branch are 
bootstrap support based on 1000 replications and values lower than 50% 


were omitted. 


3 Discussion 


Tree shrews are received more and more attention 
due to their potential use as animal models for human 
disease and as alternatives for primates in biomedical 
research (Cao et al, 2003; Peng et al, 1991). Although the 
exact phylogenetic position of tree shrews within the 
Euarchontoglires remains of debate (Arnason et al, 2002; 
Murphy et al, 2001; Nie et al, 2008; Olson et al, 2005; 
Peng et al, 1991; Sargis, 2004; Schmitz et al, 2000), 
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available evidence has undoubtedly suggested a closer 
affinity of tree shrews to well-known primates such as 
rhesus monkeys (Macaca mulatta) and crab-eating 
macaques (Macaca fascicularis) than mice or rats. Given 
that tree shrews are also small animals and easy to raise, 
it may be promising to develop them as replacements for 
monkeys and other endangered non-human primate 
species. 

There are 15 recognized species and 40 subspecies 
in genus Tupaia (http://www.bucknell.edu/msw3/; Helgen, 
2005). However, to our knowledge an insufficient 
number of molecular studies have been done and 
necessary information for biomedical research and 
breeding as an experimental animal remain unknown. 
Work on gene expression in the hippocampal tissue of 
tree shrews treated with cortisol (Alfonso et al, 2004) 
and proteome analysis of aflatoxin Bl-induced 
hepatocarcinogenesis (Li et al, 2008) has yielded useful 
information, but compared to the available knowledge 
for mice and rats there is still much unknown about these 
animals. 

Here we focused on mtDNA sequence variation in 
Chinese tree shrews captured from urban Kunming. Our 
initial aim was to better understand the genetic structure 
of this species and to establish a matrilineal marker for 
tracing lineages in a subsequent tree shrews breeding 
project. In addition, we wanted to better understand 
overall genetic diversity of this animal resource. By 
comparing the mtDNA sequence differences between 
Chinese tree shrews and northern tree shrews, we found 
that they had a large divergence (Figs.1, 2). Network 
analyses showed that Chinese tree shrews and northern 
tree shrews are clearly separated from each other (Fig. 3). 
Our results raise caution over the routine use of tree 
shrews as animal models and genetic differences are 
ignored, especially when studies simply named their 
subjects as Tupaia belangeri. 

The divergence of Chinese tree shrews and northern 
tree shrews was unexpected and implies that a systematic 
analysis of the six proposed subspecies of Chinese tree 
shrew (Wang, 1987) requires verification in future 
studies. Moreover, capturing and domesticating wild 
animals with the representative lineages should be 
encouraged during our efforts in establishing inbred tree 
shrew strains for medical and biological research. The 
observation that the tree shrew population inhabiting 
Kunming has a considerably high mtDNA genetic 
diversity is consistent with the result of a recent study of 
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microsatellite markers in which He et al (2009) found 
that nine out of 11 loci exhibited polymorphisms in 30 
Chinese tree shrews. Data such as that presented here 
will be helpful in selecting founding individuals for the 
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